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Abstract 

We have studied the scalar i;^^-model in the symmetric phase and the non- 
compact U{1) gauge theory on a momentum lattice using the Langevin equation 
for generating configurations. In the (/)^-model we have analyzed the renormal- 
ized mass and in the C/(l)-model we have analyzed the Wilson loop operator. 
We used a second order algorithm for solving the Langevin equation, and we 
looked for the convergence rate of the method. We studied the stochastic time 
needed to generate equilibrium configurations and compared first and second 
order schemes for both models. 



1 Introduction 



One major problem in lattice gauge theory has been the increased amount of computer 
time and memory needed to do simulations near a critical point. The so-called critical 
slowing down can be countered, e.g., with the Fourier acceleration method Many 
other methods have been developed to accelerate convergence of the calculation, e.g., 
the cluster algorithm p[, suitable for Ising-type systems. When we studied the scalar 
0^-model near the critical point on a momentum lattice 0, it turned out that quite 
reasonable results for the renomalized mass, wave function renormalization, etc., could 
be obtained on relatively small lattices (3^ — 7^) . In this note we want to report a more 
detailed numerical analysis of the convergence behavior on a momentum lattice. 



2 Momentum Space 

Momentum space has several interesting advantages over coordinate space: 

• The kinetic energy of the action is local. 

• One can implement Fourier acceleration to fight critical slowing down. 

• At the critical point, the correlation length goes to infinity. The behavior at 
X oo in coordinate space corresponds to — in momentum space. 

• The correlation function in momentum space behaves as l/(m|j + k"^) near the 
critical point. At the critical point, 0, but one can stay away from the 
pole by choosing some k"^ ^ 0. 

• Although the action in momentum space is non-local, one can use a fast Fourier 
transform to switch to coordinate space, where the interaction is local. 



3 Stochastic Quantization 



We have used stochastic quantization, which is a direct method to go from the contin- 
uum formulation to practical algorithms for numerical simulations. We have studied 
the model and the non-compact f/( 1)3+1 gauge model. Stochastic quantization 
was introduced by Parisi and Wu [Q, and it has been applied by many authors (for a 
review see Damgaard and Hiiffel Hi). The idea is to consider the Euclidean quan- 
tum field as the equilibrium state of a statistical system coupled to a heat reservoir. 
The evolution of this statistical system is described by the Langevin equation, which 
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reads in momentum space, 



or S(j){-k,T) 

where S[(f)] is the action in momentum space and fjlk, r) is a field of Gaussian noise. 
When r — > oo the statistical system is "equal" to the Euclidean field; r was introduced 
as fictitious time to evolve the system. Our goal is to solve the Langevin equation 
many times to generate a good representation of the canonical ensemble. Then we 
perform the canonical averaging of the desired observable. 



4 Numerical Point of View 

The expectation value of an observable is defined by 

= lim - rdT'OWr')]. 

Because we suppose that our system is ergodic, the expectation value can be computed 
over the configuration space, 

1 ^ 

n=l 

So one must solve the Langevin equation N times in parallel, which can be done using 
a parallel processor. 

The Langevin equation is a first order differential equation that we integrate up to 
a final value Tfinai, by increments of At. It can be solved using a standard method like 
Euler's method where the truncation term is O(Ar^), or a higher order scheme like 
Heun's method (second order Runge-Kutta) where the truncation term is O(Ar^). 
Batrouni et al. IQ have shown that a higher order integration scheme has many 
advantages, and our experience agrees with theirs. For the same level of error, 
Heun's method is faster; it takes fewer "time" steps. For the same "time" step, it was 
46% slower for the 0'^-model, but in the case of the non-compact f/(l)-model it was 
only 12% slower; in the case of the 0^-model, the second order scheme required the 
evaluation of two additional fast Fourier transforms, but in the case of the f/(l)-model 
only a few more additions and multiplications were required. In Fig. [1] we display for 
the 0^-model the behavior of the renomalized mass m^j in the neighborhood of the 
critical point. We have computed in Fig. [2] the renormalized mass as a function of 
Tfinai for large At on a lattice of 3"^ points. This corresponds to a coupling parameter 
K = 0.12025, while the critical point is at Kc = 0.1257(1). Stability of the solution is 
obtained for a value of Tf^ai ~ 12.5 by Heun's method. For larger values of Tfinai, the 
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accumulation of error becomes important and fluctuations can appear. We can see 
that the value of itir drops beyond Tfmai ~ 25 in Fig. [2]. It is easier to determine the 
stability region in the second order method. 

For the non-compact ?7(l)-model we have evaluated the expectation value 
Yll,,u=iFfiuFtJ.u'), where F^^, is the field tensor. We get very good agreement with 
the known exact value; an error of less than 1% can be obtained easily. In momentum 
space {F^yF^y) is a local observable. 

We also computed the expectation value W{I , J) of the Wilson loop observable, 
which corresponds to the creation of a quark- ant iquark pair, the propagation of this 
"meson" and the destruction of the pair (/ x J being the size of the rectangular 
loop measured in units of the lattice spacing). It is a non-local observable that can be 
computed analytically in the f/(l)-model. It seems to be more sensitive to fluctuations 
than {F^^F^^). When integrating the Langevin equation, if we continue to iterate too 
long after the equilibrium is reached, the field of noise fjik, r) creates some fiuctations 
in the solution. These fluctuations are more important for non-local observables, 
which are difficult to evaluate on small lattices. In Figs. [3] and [4] , we computed the 
expectation value of the Wilson loop operator for different sizes of loops and for two 
different "time" steps on a lattice of 6^ points. We can see the fluctuation introduced 
by the noise field T){k,T). In Fig. [3] the "time" step was Ar = r/j„a//100, and in 
Fig. [4] the "time" step was Ar = r/j„„//150. The important difference is the number 
of "time" steps. When we compare the results of Figs. [3,4] we notice a difference 
which seems to exceed the error (9(Ar^) due to finite Ar. We can reduce fiuctuations 
by increasing the number of configurations of fields, but we can also average several 
values of r up to t final ■ Although we have not used the last method, we think it might 
reduce fiuctuations 0. 

Stochastic quantization allows us to choose the value of Tfi^ai- The propagator for 
the free scalar theory in Euclidean momentum space for finite Tfinai is 

(mm) = S{k + ,>^ ^-^M-2rf.Uf + m^)) ^ 
\ / fc^ -|- m"^ 

compared with the Euclidean Feynman propagator, 

(^{k^{k')) = S{k + k')—^^. 

We can see that one must choose ^Tjinaiik? + m^)) large enough to make the expo- 
nential term negligible. For the non-compact f/(l)-model in momentum space, we are 
free to choose a value for each lattice point (Fourier acceleration). The point /c = 
turns out to be dangerous in the ?7(l)-model because there is no interaction between 
lattice sites. For the 0'^-model that does not seem to be the case; the convergence 
rate seems to be a little slower near the critical point, but does not differ very much. 
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5 Conclusion 



Our study of two models shows the feasibihty of doing lattice simulations on a small 
momentum lattice. The second order scheme was by far more appropriate because of 
its greater accuracy and speed. The difference in the results due to variation of At 
encountered during the simulation of the f/(l)-model for the Wilson loop operator 
does not seem to pose a major problem. We believe that the averaging over several 
values of r up to Tfinai can increase the numerical stability. This kind of behavior has 
not been seen in the 0^-model. The knowledge of the stability region is very important 
for simulations of theories like SU{2) and SU{3). For these models, the time needed 
for generating equilibrium configurations is much larger than for the C/(l)-model. 
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Figure Captions 

Fig.l Renormalized mass rriR of the 0^-model versus coupling parameter k. The 
semi-analytical results are taken from Liischer and Weisz Ref. [0. The full 
curve is a fit to Liischer and Weisz's data to guide the eye. 

Fig. 2 Renormalized mass mn of the 0^-model versus stochastic time Tfinai- 

Fig. 3 Wilson loop W{I, J) of the t/(l)-model versus coupling constant g; At = 

Tfinai/ 100. 

Fig.4 Same as Fig. [3] but At = Tfinai/150. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/9310012vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/9310012vl 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/9310012vl 



This figure "figl-4.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/9310012vl 



